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applications 

Silvana Botti,^'^ Jose A. Flores-Livas,^ Maximilian Amsler,^ Stefan Goedecker,^ and Miguel A.L. Marques^ 

^ Laboratoire des Solides Irradies and ETSF, Ecole Poly technique, CNRS, CEA-DSM, 91128 Palaiseau, France 

^ Universite de Lyon, F-69000 Lyon, France and LPMCN, CNRS, UMR 5586, Universite Lyon 1, F-69622 Villeurbanne, France 

^Department of Physics, Universitdt Basel, Klingelbergstr. 82, 4056 Basel, Switzerland 

(Dated: March 1, 2013) 

We present state-of-the-art first-principle calculations of the electronic and optical properties of 
silicon allotropes with interesting characteristics for applications in thin- film solar cells. These 
new phases consist of distorted sp"^ silicon networks and have a lower formation energy than other 
experimentally produced silicon phases. Some of these structures turned out to have quasi-direct and 
dipole-allowed band gaps in the range 0.8-1.5 eV, and to display absorption coefficients comparable 
with those of chalcopyrites used in thin- film record solar cells. 
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If Nature chose carbon as the scaffolding of life, 
mankind chose silicon as the building block of much of 
the high-technology of today. In fact, its advantages are 
many: silicon is the second most abundant element in 
the Earth's crust, its processing is well controlled at the 
industrial level, and its electronic properties and defect 
physics are theoretically well understood. As an elemen- 
tal substance, it is an intrinsic semiconductor that can 
readily be p- and n-doped with a multitude of different 
elements. This makes silicon the material of choice for 
applications in electronics. 

Silicon is also the leading player in the important field 
of photovoltaic energy production. In fact, the large ma- 
jority of the first-generation solar cells, based on a single 
p-n junction, use bulk silicon as absorber layer. Its band- 
gap {Eg = 1.12 eV at room temperature^) lies in the opti- 
mal interval of values to maximize energy conversion effi- 
ciency, according to the Schockley-Queisser limit. ^ How- 
ever, the fact that the gap is indirect, and that the op- 
tical gap is larger than 3eV, makes silicon a very bad 
absorber of sunlight^ Therefore, the absorbing layer has 
to be at least 100 /im thick and, as a consequence, the 
crystal needs to be very pure, so that the mean- free path 
of the carriers is comparable with the size of the layer 4a^ 
This is a limiting factor in the reduction of the cost of 
photovoltaic modules. 

Due to these well-known limitations, materials with 
better absorption coefficients in the visible have been 
put forward in the past years in order to replace silicon 
in thin-film solar cells. ^^ Emphasis has obviously been 
given to direct band-gap materials, with gaps between 
1 and 1.5 eV and with absorption spectra that strongly 
overlap with the solar spectrum. ^^^^ However, and in 
spite of many spectacular advances, none of these ma- 
terials has been able to dethrone silicon. Of course, the 
ideal solution would be to engineer silicon such that it 
absorbs strongly in the visible, in order to use thin- film 
technology and at the same time maintain the silicon- 
based processes that are currently employed. This would 
allow for thinner, flexible, and cheaper silicon solar cells. 

Aiming at maximizing the spectral overlap, some au- 
thors have recently suggested different ways to manipu- 



late the optical properties of silicon. Nanostructuring is 
a largely explored way to obtain direct-gap silicon by re- 
laxing translational symmetry J^ However, the existence 
of a direct gap does not guarantee that transitions at the 
absorption edge are dipole allowed. Indeed, this is usu- 
ally not the case in silicon nanostructuresJ^ D'Avezac et 
ali^ proposed recently an ultrathin silicon/germanium 
superlattice with excellent absorption properties. The 
experimental synthesis of such system requires, however, 
control of the growth of pure monolayers. Compensated 
doping of silicon obtained through substitutional impu- 
rities was also shown to increase by 25% the absorption 
of 10 /im thick silicon layers. ^^ Nevertheless, doping is 
expected to introduce recombination centers detrimental 
for the performances of the device. 

Other works investigated the possibility to modify the 
band structure of silicon by using allotropic phases with 
different crystal symmetries J^— The most stable phase 
of silicon, at ambient conditions, is the cubic diamond 
structure. High pressure phases have been experimen- 
tally studied up to 248 GPa^ and several calculations of 
the phase transitions are present in literature j^^Ti^i Upon 
increase of pressure silicon exhibits a series of phase tran- 
sitions: from cubic diamond to /3-Sn at around 12 GPa, 
from /3-Sn to orthorhombic (Imma symmetry)^ and then 
to simple hexagonal at 13-16 GPa, from simple hexagonal 
to an orthorombic Cmca phase at about 38 GPa,^^ from 
Cmca to hexagonal close pack at 42 GPa, and finally to 
face-centered cubic at 78 GPa. It has been known for 
more than 20 years that, upon slow release of pressure 
from the /3-Sn structure, silicon transforms into the rom- 
bohedral R8 and the body-centered BC8 crystals, char- 
acterized by distorted sp^ bondings j^^Ti^i R8 and BC8 
are met ast able phases, as they keep existing at ambient 
conditions. If BC8 silicon is heated to temperatures in 
the range from 200 to 600 °C, it transforms to another 
metastable phase: the lonsdaleite hexagonal diamond^ 
Two other phases, the so-called Si- VIII and Si-IX struc- 
tures, were observed experimentally upon rapid release 
of pressure from /3-Sn silicon. ^^ However, their crystal 
lattices were not fully characterized. Nanoindetation ex- 
periments gave evidences of another phase, known as Si- 



XIII^ although also in this case very little information 
is available on its crystal structure. Other met a-st able 
phases that lay close in energy to cubic diamond sili- 
con have been proposed theoretically, using a multitude 
of methods. ^1^— Malone et at. calculated the band 
structure and the absorption spectrum one of these low- 
energy arrangements, charachterized by a body-centered- 
tetragonal (BCT) unit cellJ^ A similar study was carried 
out by the same authors also for the experimental R8 
structure. ^^ Unfortunately, BCT silicon did not show en- 
hanced absorption properties in the visible if compared 
to cubic diamond silicon and R8 turned out to have a too 
small gap (0.24 eV) for photovoltaic applications. How- 
ever, the low-pressure phase diagram of silicon is still 
relatively unexplored, which makes us believe that there 
might exist unexplored silicon phases with optical prop- 
erties suitable for thin- film photovoltaics. 

In this paper we present the results of a structural 
prediction search that leads to a number of previously 
unknown, low-energy, sp^ phases of silicon that have 
excellent properties for photovoltaic applications, with 
quasi-direct band-gaps between 1-1.5 eV, and excellent 
absorption properties for solar light. Note that the band 
gap imposes an upper bound on the open-circuit voltage 
of the device, implying that too small indirect band-gap 
compounds should be discarded. All structures we se- 
lected have a total energy higher than the one of cubic 
silicon by less than 0.15 eV per atom. This restriction as- 
sures that the experimental synthesis of these structures 
is energetically plausible. 

The silicon crystalline arrangements under investiga- 
tion were found by performing a structural relaxation of 
low-enthalpy carbon primitive cells with 8 atoms, which 
we had previously calculated^ by applying a very effi- 
cient structural prediction algorithm, the minima hop- 
ping method (MHM).^^'^^ The MHM was designed to 
obtain the low-enthalpy structures of a system given 
solely its chemical composition. The enthalpy surface 
is explored by performing consecutive short molecular 
dynamics escape steps followed by local geometry relax- 
ations taking into account both atomic and cell variables. 
The initial velocities for the molecular dynamics trajecto- 
ries are chosen approximately along soft mode directions, 
thus allowing efficient escapes from local minima and 
aiming towards low energy structures. The predictive 
power of this approach has already been demonstrated 
for a wide range of applications.—!^— As carbon ad- 
mits both sp^ and sp^ bonds, its phase diagram is much 
richer than the one of silicon. For this reason the low- 
enthalpy carbon allotropes found at low pressure provide 
a good structural database for silicon analogues. Indeed, 
we observed that the carbon structures with sp^ bonds 
turn out to be unstable and relax into a different geom- 
etry with sp^ arrangement. 

The evaluation of energies and forces required for the 
MHM were performed within density functional theory 
(DFT) using the ABINIT code^ with simulation cells 
of 8 carbon atoms at 15 GPa. The lowest energy sili- 
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FIG. 1. Representation of the most important sp^ Si struc- 
tures found in this work. The superscript in Imma*^^-* serves to 
distinguish this structure from the high-pressure phase Si-XI 
of the same symmetry. 



TABLE I. Selected low-energy silicon allotropes found in this 
study. Total energies per atom are given relatively to the 
diamond structure. PBE indirect gaps are compared with 
corresponding GW indirect gaps. All energies are in eV. The 
references to works that present carbon analogues are given 
in the last column. 



structure 


spc group 


^totai/atoi 




T^GW 

^gap 


C analogue 


Z 


65 


0.06 


0.72 


1.22 


Ref. 34,^, and^ 


M 


12 


0.07 


0.46 


0.99 


Ref. 56 


M-10 


10 


0.08 


0.97 


1.40 


Ref. 45 and 57 


C222i 


20 


0.08 


0.90 


1.39 


Ref. 59 


Imma^^) 


74 


0.12 


0.63 


1.15 


- 


Cmcm 


63 


0.13 


0.32 


0.83 


Ref. 57 


P-1 


2 


0.13 


0.28 


0.75 


- 


P2i/c 


14 


0.14 


0.73 


1.29 


- 



con analogues were further relaxed and characterized us- 
ing the VASP code.^^ We selected our /c-point grids to 
ensure an accuracy of 0.01 eV in the total energy, and 
all forces were converged to better than 0.005 eV/A. To 
approximate the exchange-correlation functional of DFT 
we used the Perdew-Burke-Ernzerhof (PBE)^ general- 
ized gradient approximation. 

In an energy interval of 0.15 eV per atom (around 
1740 K) from the diamond structure, and ignoring the 
well-known cubic diamond and hexagonal diamond sili- 
con, we found a total of 16 structures. We checked the 
stability of these phases by calculating their phonon band 
structure, in order to make sure that all phonon modes 
are real. We note that the experimentally known meta- 
stable phases of silicon, R8 and BC8, have higher ener- 
gies than the phases considered here. Interestingly, also 
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FIG. 3. (Color online.) Absorption spectra of the most im- 
portant sp^ Si structures found in this work compared to the 
reference air mass 1.5 solar spectral irradiance,— given in ar- 
bitrary units. 



FIG. 2. (Color online.) Energy per atom as a function of 
volume per atom for the most important sp^ silicon structures 
found in this work. The zero of energy is the total energy 
per atom of the cubic diamond structure in equilibrium. For 
comparison we also include some already known structures of 
silicon, e.g., cubic diamond, a clathrate,— BC8, R8. 



the most studied theoretical meta-stable phases (ST12, 
BCT, Ibam etc.) are energetically less favored than the 
new structures presented in this work. 

All our structures have sp^ bonding, and are semi- 
conducting with Kohn-Sham PBE (indirect) gaps rang- 
ing from 0.3 to 1.2 eV. Some of the carbon analogues 
can already be found in the literature (see Table H]). It 
is noticeable that none of the around 300 semiconduct- 
ing structures identified by our simulation^^ has a di- 
rect quasi-particle gap. Nevertheless, several of them ex- 
hibit a quasi-direct gap. The lowest-energy structures 
we found with promising quasiparticle gaps are listed, 
together with their space group number and the value of 
the calculated gaps, in Table [H The corresponding atomic 
arrangements are depicted in Fig. [TJ 

In Fig. [2] we show the total energy per atom as a func- 
tion of volume per atom for the new silicon allotropes 
found in this study. For the sake of comparison we in- 
cluded the energy curves of some well known phases of 
silicon, e.g., cubic diamond silicon, a clathrate, the BC8 
and R8 structures. 

By inspection of Fig. [2] and Table |I] we can easily de- 
duce that, among the novel silicon structures, the lowest 
energy allotropes is Z-silicon, followed by M and M-10 
silicon. These crystal structures have carbon counter- 
parts that are thought to play an important role in cold- 
compressed graphite. ^^ The volumes per atom at equi- 
librium of the new structures are between 20 and 22 A^, 
which places them in an intermediate position between 



the open clathrates and the experimentally observed BC8 
and R8 structures. 

In order to characterize the fundamental and optical 
band gaps, together with the absorption spectra, we per- 
formed perturbative GW— calculations and we solved 
the Bethe-Salpeter equation^i^ using the abinit pack- 
aged We want to stress that we employed the most 
accurate methods available in the community to study 
electronic excitations, that have proved to give excel- 
lent results for a wide range of sp semiconductors and 
insulators^i^ In particular, the BSE absorption spec- 
trum of cubic diamond silicon— is in excellent agreement 
with experimental measurements^^ thanks to the inclu- 
sion of excitonic effects. We used Troullier-Martins norm- 
conserving pseudopotential^ and the PBE exchange- 
correlation functional!^ The planewave cut-off energy 
for all runs was 15Hartree. For GW calculations the 
cutoff of the dielectric matrix was set to 5 Hatree and we 
needed 12 states per atom, together with the method of 
Bruneval and Gonze, to achieve a convergence of around 
0.05 eV in the band gap. For the Bethe-Salpeter calcu- 
lations we used a cutoff of 2.5Hartree for the dielectric 
matrix and a k-point spacing of around 0.012 in recipro- 
cal lattice vector units for the shifted k-point mesh. This 
criterion yields a 14x14x14 shifted grid for the diamond 
structure. 

All the GW corrections to the Kohn-Sham PBE band 
structures can be approximated by a rigid shift of the 
conduction bands, as they are fairly independent of the 
k-point s. Furthermore, this rigid shift does not change 
significantly among the different crystal lattices, ranging 
from 0.4 to 0.6 eV, even if the Kohn-Sham gaps exhibit 
larger variations. By inspecting the values for the in- 
direct gaps in Table |I] we can observe that several sili- 
con structures have a gap in the optimal frequency in- 
terval for maximizing photovoltaic efficiency. However, 
one should not forget that all the gaps that we calcu- 



lated are indirect, even if in many cases the direct gaps 
are much smaher than in the cubic diamond phase. This 
fact suggests that it is conceivable that the absorption 
edge for direct transitions of these metastable ahotropes 
is at significantly lower energies than in conventional sili- 
con. However, information on the band structure alone is 
not conclusive, since the dipole matrix elements between 
states close to the valence band maximum and the con- 
duction band minimum can be very small and suppress 
light absorption close to the absorption edge. In view of 
that, we performed accurate calculations of the absorp- 
tion spectra by solving the Bethe-Salpeter equation. 

In Fig. [3] we show the calculated absorption spectra 
of the allotropes that are most promising for applica- 
tions in solar cells. The calculated absorption spectrum 
of the cubic diamond phase is also shown for compari- 
son, together with the reference air mass (AM) 1.5 solar 
spectral irradiance. Note that indirect absorption con- 
tributions, which arises from phonon assisted interband 
transitions are not included in our calculations. All new 
structures start to absorb very close to their direct gap, 
with the result that their absorption spectra overlap sig- 
nificantly with the solar spectrum. Remarkably, the ab- 
solute optical absorption between 1.5 and 3eV strictly 
is comparable to that of Cu(In,Ga)Se2 compounds, that 
are regarded as excellent absorbers for thin-film solar cell 
technology. ^^ This indicates that these new silicon phases 
could potentially be employed in highly efficient silicon 
thin-film solar panels. 

If we make a critical evaluation of all calculated proper- 



ties, we can conclude that the most promising allotropes 
are the M-10 and the C222i. In fact, their total ener- 
gies per atom are particularly low, their indirect gap is 
1.4 eV, while the direct gap is only slightly larger (1.5 eV 
for M-10 and 2.0 eV for C222i). Moreover, absorption is 
dipole allowed starting from the absorption edge, as it is 
revealed by the long tails of their BSE absorption spectra 
in the visible. 

In conclusion, we predict several novel metastable 
phases of silicon, obtained through a structural search 
based on carbon analogues found using the minima hop- 
ping method. We filtered the most promising structures 
for use in thin-film photovoltaics by imposing that the to- 
tal energy per atom does not differ by more than 0.15 eV 
from the total energy per atom of cubic diamond silicon. 
Moreover, we imposed that the fundamental band gap is 
in the range 1.0-1.5 eV. Some of the new structures re- 
vealed a strong absorption in the visible, with absorption 
coefficients comparable to those of chalcopyrite absorbers 
used in thin-film record solar cells. These results call for 
further experimental studies of the low-pressure phase 
diagram of silicon, and could open new routes to design 
highly efficient thin-film silicon solar cells. 
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